Effects of spin-orbit interaction on the magnetic and electronic structure of 

antiferromagnetic LaFeAsO 
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Magnetic and electronic structures in LaFeAsO in the single-stripe-type antiferromagnetic (AFM) 
phase are studied using first-principles density-functional calculations including the spin-orbit in- 
teraction. We find that the longitudinal ordering (LO) where Fe magnetic moments are parallel 
or anti-parallel with the in-plane AFM ordering vector, q, is lower in energy than transverse or- 
derings (TOs) where Fe magnetic moments are perpendicular to q, in excellent agreement with 
recent neutron diffraction experiments. Calculated energy difference between LO and TOs is about 
0.2 meV per Fe atom, indicating that LO will prevail at temperature below 2 K. We also find that 
the spin-orbit interaction splits degenerate bands at k = (l/2)q by about 60 meV, depending on 
spatial directions of the Fe magnetic moments. 

PACS numbers: 74.70.Xa, 71.15.Mb, 75.70.Tj, 71.70.Ej 



Superconductivity and magnetism in iron oxypnictides 
LaFeAsO and other i?FeAsO (R = rare-earth) have be- 
come central issues since their discovery* 1 - - — These ma- 
terials exhibit superconductivity at high temperatures, 
for example, 55 K in SmFeAsOi-^F^ji Without dop- 
ing, these compounds show antiferromagnetism at low 
temperature, and they become superconductors when 
doped with electrons or holes. The occurrence of su- 
perconductivity in the vicinity of magnetic phase sug- 
gests that understanding the magnetic properties could 
be essential in revealing the origin of their unconventional 
superconductivity^ — 

As temperature lowers, undoped LaFeAsO undergoes a 
structural transition from a tetragonal to an orthorhom- 
bic structure and Fe magnetic moments form a single- 
stripe-type antiferromagnetic (AFM) ordering, where Fe 
magnetic moments are ordered antiferromagnetically in 
an orthorhombic in-plane axis, say a-axis, ferromagnet- 
ically in the other in-plane axis, say o-axis, and antifer- 
romagnetically along the c-axis^ Density functional cal- 
culations have been successful in predicting and explain- 
ing the in-plane AFM orderings of Fe magnetic moments 
in iron pnict ides^^ i 10 i 12 i 15 i 16 and chalcogenide o 17 ' 18 cor- 
rectly although the calculations predict substantially 
large values of Fe magnetic moments. 

More recently, neutron diffraction experiments 1 ^ on 
the magnetic structure in LaFeAsO claimed that spatial 
directions of Fe magnetic moments in the single-stripe- 
type AFM phase are collinear with the AFM direction in 
the Fe plane, i.e., the orthorhombic a-axis, rather than 
other directions^ This result suggests a significant role 
of the spin-orbit interaction in LaFeAsO, since magnetic 
moments are coupled to the lattice structure via the spin- 
orbit interaction only; yet, there is no theoretical report 
on this issue. 

In this Communication, we will present electronic 
and magnetic properties in LaFeAsO calculated by first- 



principles density-functional methods including the spin- 
orbit interaction. We will show that spatial directions of 
Fe magnetic moments in the ground state of the single- 
stripe-type AFM are along the in-plane AFM direction, 
i.e., the orthorhombic a-axis, which is consistent with 
the recent neutron diffraction experiments^ We will also 
show that the energy cost for aligning Fe magnetic mo- 
ments perpendicular to the orthorhombic a-axis is about 
0.2 meV/Fe, indicating that the Fe magnetic moments 
prefer the orthorhombic a-axis direction significantly at 
temperatures below 2 K. In addition, we show that the 
spin-orbit interaction splits degenerate bands at a high- 
symmetry point, X, in the orthorhombic Brillouin zone 
by about 60 meV, depending on spatial directions of the 
Fe magnetic moments. 

Our first-principles calculations are performed with 
the generalized gradient approximation 2 ^ to the den- 
sity functional theory and the projector-augmented-wave 
pseudopotential o 21 i 22 as implemented in VASP i 23 ' 24 We 
regard 4p 6 5s 2 6d 1 electrons in La, 3p 6 3d 6 4s 2 in Fe, 4s 2 4p 3 
in As, and 2s 2 2p 4 in O as valence electrons. All of our cal- 
culations are performed using a 8 x 8 x 4 Monkhorst-Pack 
/c-point mesh2£ centered at T for orthorhombic (Cmme) 
structure. Electronic wavefunctions are expanded with 
plane waves up to a kinetic-energy cutoff of 500 eV, and 
the tetrahedron method with the Blochl corrections is 
used for the Brillouin zone integration. 2 - With the spin- 
orbit interaction included in the total-energy functional, 
we carefully check the convergence of the self-consistent 
calculation of the electron density in order to distinguish 
small difference in the total energy caused by difference 
in the spatial directions of the Fe magnetic moments. 

We optimize atomic positions in LaFeAsO using a 32- 
atom orthorhombic supercell, a x b x 2c, which is twice 
the structure shown in Fig. 1(a) along the c-axis. Dur- 
ing the optimization, the supercell lattice constants are 
fixed to a = 5.7063 A, b = 5.6788 A, and c = 8.7094 A, 
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(a) LaFeAsO (b) LO 




FIG. 1: (Color online) Atomic and magnetic structures in 
LaFeAsO. (a) An orthorhombic unit cell of LaFeAsO in the 
single-stripe-type antiferromagnetic phase. The unit cell con- 
tains four formula units, i.e., 16 atoms, (b) Longitudinal or- 
dering (LO), (c) in- plane transverse ordering (TO), and (d) 
out-of-plane TO of Fe magnetic moments in the FeAs layer. 
In LO, in-plane TO, and out-of-plane TO, spatial directions 
of Fe magnetic moments are parallel or anti-parallel with the 
orthorhombic a, b, and c directions, respectively. In (b) and 
(c), directions of Fe magnetic moments are denoted by arrows. 
In (d), Fe magnetic moments in the upper and the lower rows 
(in red) are in the c direction, while those in the middle row 
(in blue) are in the opposite direction. 



which are the measured values at 2 K,— and the single- 
stripe-type AFM order is imposed, where Fe magnetic 
moments are antiferromagnetic along the a direction and 
ferromagnetic along the b direction. Atomic positions are 
fully optimized until residual forces on them are less than 
0.005 eV/A. The spin-orbit interaction is considered dur- 
ing the structural optimization, but obtained equilibrium 
positions of atoms depend negligibly on spatial directions 
of Fe magnetic moments as long as they are in the single- 
stripe-type AFM. 

To investigate effects of the spin-orbit interaction on 
the Fe magnetic-moment directions, we consider three 
different configurations, as shown in Figs. l(b)-(d), where 
the Fe magnetic moments are along the orthorhombic a, 
b, and c directions, respectively. In Fig. 1(b), the Fe 
magnetic moments are in a longitudinal order (LO) in 
the sense that they are parallel or anti-parallel to the 
in-plane AFM direction, i.e., the a-axis in our supercell. 
In Figs. 1(c) and (d), the Fe magnetic moments are in 



TABLE I: Calculated total energies, E to t, and Fe magnetic 
moments mF Cl in LaFeAsO with various spin-configurations, 
obtained including the spin-orbit interaction. Considered spa- 
tial directions of Fe magnetic moments are the longitudinal 
ordering (LO), the in-plane transverse ordering (TO), and the 
out-of-plane TO, as shown in Fig. 1. For each spin configura- 
tion, we consider antiferromagnetic (AFM) and ferromagnetic 
(FM) orderings along the c axis; in the former Fe moment di- 
rections are in phase along the c axis, while in the latter they 
are out of phase along the c axis. The total energy of LO with 
AFM along the c-axis is set to zero. Without the spin-orbit 
interaction, all spin-configurations have the same total energy 
and the same Fe magnetic moment, 1.689 /^b/Fb. 





spin 


-Etot 


TTlFc 




configuration 


(meV/Fe) 


(Ms/Fe) 


AFM 


LO 


0.00 


1.680 


along 


in- plane TO 


0.25 


1.682 


c-axis 


out-of-plane TO 


0.15 


1.681 


FM 


LO 


0.00 


1.680 


along 


in- plane TO 


0.25 


1.682 


c-axis 


out-of-plane TO 


0.15 


1.681 



in-plane and out-of-plane transverse orders (TOs), re- 
spectively, in the sense that they are perpendicular to 
the in-plane AFM direction and they are either in the 
Fe plane or out of the plane. Although the spatial direc- 
tions of Fe magnetic moments are different in the three 
configurations, all of them have the same single-stripe- 
type AFM in the sense that neighboring Fe magnetic 
moments are in opposite directions along the a direc- 
tion while they are in the same direction along the b 
direction. For each of the three configurations, we con- 
sider two different inter-Fe-plane magnetic orderings: fer- 
romagnetic and antiferromagnetic orderings between Fe 
magnetic moments neighboring along the c direction. 

In our density-functional calculations, we do not im- 
pose any constraint on the spatial direction of the Fe 
magnetic moments; however, when we take one of the 
three configurations as initial guess for the spin density 
at the start of the self-consistent iteration in our cal- 
culations, the spatial directions of the Fe magnetic mo- 
ments are stationary during the iteration until the self- 
consistency is reached. In the cases that we try spatial 
directions of Fe magnetic moments slightly away from the 
LO case as initial guess, the directions of the Fe magnetic 
moments are changed during the self-consistent iteration 
and they converge to the LO case. 

Table I shows our calculational results for the total 
energy and Fe magnet moments of LaFeAsO including 
the spin-orbit interaction. We find that LO is the lowest- 
energy configuration, while in-plane TO and out-of-plane 
TO are higher in energy than LO by 0.25 meV and 
0.15 meV per Fe atom, respectively. The total energies in 
our calculations are independent of c-axis orderings of Fe 
magnetic moments because of rather large distance be- 
tween Fe planes in LaFeAsO. Obtained Fe magnetic mo- 
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FIG. 2: (Color online) Calculated spin densities in LaFeAsO 
in the single-stripe-type AFM phase, represented with iso- 
surfaces. (a) The a-axis component of the spin density with 
the spin-orbit interaction in the case of the longitudinal order- 
ing (LO) of Fe magnetic moments, (b) Spin density without 
the spin-orbit interaction, drawn for comparison. In (a) and 
(b), left panels show top views of the spin density in Fe planes 
in which Fe atoms are at corners, at middles of edges, and at 
the center of the square, and the right panels show corre- 
sponding three-dimensional views of the spin density in the 
orthorhombic unit cell. In each top view, iso-surfaces in the 
upper and lower rows (in red) correspond to a positive spin 
density while those in the middle row (in blue) correspond to 
a negative one. 



merits are 1.680 /Lte/Fe for LO, 1.682 hb/Fc for in-plane 
TO, and 1.681 /^e/Fe for out-of-plane TO, which are also 
independent of c-axis orderings of Fe magnetic moments. 
Without the spin-orbit interaction, all the three config- 
urations have the same total energy and the same Fe 
magnetic moment, 1.689 /ie/Fe. 

Our result that LO is the lowest-energy configuration 
of Fe magnetic moments in LaFeAsO is in excellent agree- 
ment with the recent neutron diffraction experiments 19 
where spatial directions of Fe magnetic moments at 2 K 
in LaFeAsO were claimed to be along the AFM direction 
in the Fe plane. The same neutron diffraction experi- 
ments report that the size of the Fe magnetic moment in 
LaFeAsO is 0.63 /ie/Fe at 2 so our calculated sizes 
of Fe magnetic moments are about one Bohr magneton 
larger than the measured value. Since our calculated en- 
ergy gain of LO is about 0.2 meV/Fe when compared 
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FIG. 3: (Color online) Electronic energy bands along the 
high symmetry lines in LaFeAsO for three different spin- 
configurations: (a) LO, (b) in-plane TO, and (c) out-of-plane 
TO. Dark gray (red) and light gray (green) lines represent 
results obtained with and without the spin-orbit interaction, 
respectively, (d) One eighth of the first Brillouin zone, where 
high symmetry points and lines are marked. In (a) and (c), 
splitting of degenerate bands at the X point by the spin-orbit 
interaction is about 60 meV. In (d), F-X, F-Y, and T-Z lines 
are along the orthorhombic a-, b-, and c-axis in Fig. 1, re- 
spectively. 



with the other magnetic configurations and this energy 
gain corresponds to a temperature of about 2 K, our re- 
sult implies that LO is thermodynamically preferred at 
temperatures below 2 K. 

As shown in Table I, obtained magnitudes of Fe mag- 
netic moments are very weakly dependent on the spatial 
directions of the moments, and the spin-orbit interaction 
produces only 0.008 /ue/Fe difference. For more detailed 
comparison, we obtain the spatial distributions of the 
spin densities in LaFeAsO with and without the spin- 
orbit interaction, as shown in Fig. 2, and find that the 
maximal difference in the spin densities in the two cases 
is about 10% of the maximal size of the spin density in 
each case. In Fig. 2, the spin densities with and without 
the spin-orbit interaction look quite similar to each other, 
but it can be seen more easily in the top views in Fig. 2 
that the shapes of the iso-surfaces are slightly different. 

We study effects of the spin-orbit interaction on the 
electronic band structures in LaFeAsO. Figure 3 shows 
the band structures in orthorhombic LaFeAsO, where 
dark gray (red) and light gray (green) lines show the band 
structures with and without the spin-orbit interaction, 
respectively. Since most of the band structures without 
the spin-orbit interaction coincide with those with the 
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spin-orbit interaction and the former is drawn prior to 
the latter, the band structures without the spin-orbit in- 
teraction, plotted in light gray (green), are visible only 
at fc-points where the spin-orbit interaction produces sig- 
nificant difference. Band structures without the spin- 
orbit interaction are consistent with previous theoretical 
results.-^ Our results show that the spin-orbit interac- 
tion splits degenerate bands at the X point in the LO case 
[Fig. 3(a)] and in the out-of-plane TO case [Fig. 3(c)], and 
the splitting is about 60 meV. However, band structures 
in the in-plane TO case [Fig. 3(b)] do not show significant 
splitting at the X point even with the spin-orbit interac- 
tion. This dependence of band-energy splitting on the 
spatial direction of Fe magnetic moments is due to spe- 
cific orbital characters of the electronic states at the X 
point. 

In summary, we have investigated the effects of the 
spin-orbit interaction in the magnetic- and electronic 
structures of the FeAs-based compound LaFeAsO using 



first-principles density functional calculations. In our re- 
sults, the longitudinal ordering of Fe magnetic moments, 
in which the moments are parallel or anti-parallel with 
the in-plane AFM direction, is the lowest-energy mag- 
netic structure among various spin-configurations. This 
is in good agreement with the experimental results^ al- 
though the magnitude of our calculated magnetic mo- 
ment is about one Bohr magneton larger than the ex- 
perimental value. We also find that the energy gain 
of the directional alignment of Fe magnetic moments in 
LaFeAsO is about 0.2 meV/Fe, equivalent to a temper- 
ature of 2 K, and that the spin-orbit interaction splits 
degenerate bands at the X point by about 60 meV, de- 
pending on the spatial direction of the Fe magnetic mo- 
ments. 
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